FNAL-PUB-93/039-A 
May 1993 



On Identifying the Present-day Vacuum Energy 
with the Potential Driving Inflation 



Richard A. Frewin 1 

Astronomy Unit, School of Mathematical Sciences, 
Queen Mary and Westfield, Mile End Road, London El 4NS, U. K. 

James E. Lidsey 2 

NASA/Fermilab Astrophysics Center, 
Fermi National Accelerator Laboratory, Batavia, IL 60510, U. S. A., and 
Astronomy Unit, School of Mathematical Sciences, 
Queen Mary and Westfield, Mile End Road, London El 4NS, U. K. 

Abstract 

There exists a growing body of observational evidence supporting a non- 
vanishing cosmological constant at the present epoch. We examine the possi- 
bility that such a term may arise directly from the potential energy which drove 
an inflationary expansion of the very early universe. To avoid arbitrary alter- 
ations in the shape of this potential at various epochs it is necessary to intro- 
duce a time-dependent viscosity into the system. The evolution of the effective 
Planck mass in scalar-tensor theories is a natural candidate for such an effect. 
In these models there are observational constraints arising from anisotropies in 
the cosmic microwave background, large-scale galactic structure, observations 
of the primordial Helium abundance and solar system tests of general relativity. 
Decaying power law and exponential potentials are considered, but for these 
models it is very difficult to simultaneously satisfy all of the limits. This may 
have implications for the joint evolution of the gravitational and cosmological 
constants. 
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1 Introduction 



The solution to vacuum general relativity with a cosmological constant A is de Sitter 
space and this constant and solution have often been invoked to reconcile theory with 
observation. Originally Einstein believed the universe to be static and introduced a 
constant A-term into his field equations to cancel the expansionary behaviour found 
when A = 0. 1 Some decades later the steady state scenario based on de Sitter space 
was developed because observations of Hubble's constant suggested the Earth was 
older than the universe itself. 2 More recently the inflationary scenario has been pro- 
posed to solve some of the problems of the hot, big bang model. 3 During inflation the 
potential energy of a quantum scalar field dominates the energy-momentum tensor 
and behaves as a cosmological constant for a finite time. 

Realistic inflationary models predict that the current value of the density param- 
eter, Q , should be very close to unity. 4 There are a number of problems associated 
with the (A = 0, Q = 1) universe which can be resolved if A ^ 0. Firstly, its age 
is to ~ 6.52/i -1 Gyr, where h is the current expansion rate in units of 100 km sec -1 
Mpc" 1 . This is very close to the age of the oldest globular clusters in the galaxy, 
tec = (13 — 15) ± 3 Gyr, if h > 0.6 as suggested by some observations. 5 If A ^ 0, 
however, the expansion rate is increased and to m &y exceed tec if ^vacuum ~ 0.8. 6 
Moreover, most dynamical determinations of O suggest f2 = 0.2 ±0.1 up to scales 30 
Mpc and the apparent inconsistency with inflation is again resolved if O vacilum ~ 0.8. 7 
Finally, the introduction of vacuum energy into the standard cold dark matter model 
of galaxy formation accounts for the extra large-scale clustering observed in galaxy 
surveys. 8 

Hence, there are a number of reasons for supposing that a cosmological constant 
may be influential at the present epoch. 9 This work investigates whether the potential 
energy that drove the inflationary expansion could be such a term. This has been 
investigated previously within the context of general relativity by Peebles and Ratra, 10 
but their models required the form of the potential to change drastically at various 
epochs and therefore suffered from an element of ad-hoc fine-tuning. If the potential 
is relevant at the current epoch, it must either have a minimum at V ^ or contain a 
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non-vanishing decaying tail. Although the first possibility is not ruled out, it requires 
severe fine-tuning, so we focus on the second. This implies that thermalization of 
the false vacuum will not proceed via rapid oscillations of the scalar field about some 
global minimum and the form of the potential must change at various epochs. In 
general relativity the potential must be sufficiently flat for the strong energy condition 
to be initially violated, but must then become steep enough for reheating to proceed. 
But the energy density of the field must redshift at a slower rate than the ordinary 
matter components at late times if the vacuum energy is to once more dominate the 
dynamics. 54 

Instead of altering the shape of the potential we extend the gravitational sector of 
the theory beyond general relativity and investigate whether inflation was driven by 
potentials which are (a) too steep to lead to inflation in general relativity and (b) do 
not contain a global minimum. We shall refer to these as steep potentials. A number 
of unified field theories lead to scalar field potentials which exhibit both of these 
characteristics. The mechanism leading to inflation is very simple. In pure Einstein 
gravity containing a single, minimally coupled scalar inflaton field, a, the strong 
energy condition is violated if the condition a 2 < V holds, i.e. the potential energy 
dominates over the field's kinetic energy. Clearly this condition must break down 
at some point as steeper potentials are considered. But a finite interval of inflation 
is possible with steep potentials if one introduces a viscosity into the inflaton field 
equation which decays as the universe expands. This viscosity slows the field down 
and can lead to inflation. As the viscosity becomes weaker, however, the inflaton's 
kinetic energy increases significantly and a natural exit from inflation proceeds as the 
expansion becomes subluminal. We identify the dilaton field which arises in scalar- 
tensor theories as the natural source of this viscosity. 

This paper is organised as follows. We survey theories that lead to inflation with 
steep potentials in section 2. In section 3 we derive expressions for the amplitudes of 
the primordial fluctuation spectra and discuss the most stringent observational con- 
straints which any viable model of this type must satisfy. These limits arise from the 
observations of large-scale galactic structure, 11 anisotropies in the cosmic microwave 
background radiation (CMBR), 12 nucleosynthesis calculations 14 ' 15 and time-delay ex- 
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periments in the solar system. Numerical results for both decaying power law and 
exponential potentials are presented in section 4 and we conclude that successful infla- 
tion based on this mechanism is unlikely for the examples considered. Some possible 
implications of this conclusion are discussed in section 5. Unless otherwise stated, 
units are chosen such that c = H = 1, and the present day value of the Planck mass 
is normalized to m P1 = 1. 

2 Inflation with steep potentials 
2.1 Scalar- Tensor theories 

A suitable source of the viscosity is the dilaton field which arises in the Bergmann- 
Wagoner class of theories. 17 This field plays the role of a time-dependent gravitational 
constant. The field equations for these scalar-tensor theories are derived by varying 
the action 



where g = detg^ u , 1Z is the Ricci scalar, h{jp) is some arbitrary function of the 
dilaton ip, U{jp) is the dilaton self-interaction and £ ma tter is the Lagrangian of a 
perfect baryotropic fluid, which we assume to be relativistic matter with an effective 
equation of state p r = p r /3. The dilaton and inflaton field equations are coupled 
and the extra viscosity on a is due to the dynamical evolution of the effective Planck 
mass in the theory. The strength of gravity is determined by the magnitude of h(ip) 
and a fraction of the inflaton's potential energy is converted into the dilaton's kinetic 
energy rather than contributing to a 2 . 

It is well known that theories of this type appear as the low energy limits to 
a number of unified field theories and can be expressed in the Jordan-Brans-Dicke 
(JBD) form with a variable tu-parameter by defining a new scalar field \I/ = h(ip). 18 
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The action ( |2.1|) becomes 





) 



matter 



(2.2) 



where 



(2.3) 



2(d/i/W 



The theory of general relativity is recovered whenever a turning point exists in the 
functional form of h(tp) (see section 3.3), because ufi?) tends to infinity and the 
dilaton's kinetic energy decouples. Although the theories ( |2.1| ) and ( |2.2j ) are equiva- 
lent mathematically, there is a philosophical difference which arises in deciding which 
function h(ip) or should be treated as the fundamental quantity. 19 In general, if 
we consider a simple form for h{ip) such as a truncated Taylor series, this leads to a 
very complicated form for uj{^>) and vice- versa. 

For example we can expand h (ip) as some power series 



for arbitrary constant coefficients a,. To lowest order h(ip) w a 0) but this corresponds 
to a 'constant' gravitational constant and is not interesting. Moreover the linear 
term may always be eliminated by a simple field redefinition, so the lowest order of 
interest is the quadratic contribution and this is simply the standard JBD theory 
with constant uj(ty). In principle, given suitable initial conditions, inflation will then 
occur as a slowly rolls down its potential causing ip to increase. Eventually higher- 
order (i > 2) terms in the expansion (|2.4j) will become important and for appropriate 
choices of «j, such as a 2 > and a% < 0, one can easily arrange for a local maximum 
in h(%p) to exist at some value ip = t/>q- At this point the theory will be identical 
to general relativity if we normalise h(ipo) = 1. It is clear that as ip approaches 
inflation will end because the gravitational friction weakens and the inflaton speeds 
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(2.4) 



i=0 
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up. This model is a chaotic version of the hyperextended scenario and proceeds via 
a second-order phase transition. 18 

Another possibility is the induced theory of gravity where h(tp) is quadratic but 
the dilaton potential is non- vanishing and contains a global minimum at ipQ. Inflation 
driven by o could occur if ijj is initially displaced from ip , but will clearly end as 
Einstein gravity is recovered and the dilaton settles into this minimum. After spon- 
taneous compactification to four- dimensions, some Kaluza-Klein theories have this 
structure, where the dilaton is identified as the logarithm of the radius of the internal 
space. 20 If monopole and Casimir effects due to non-trivial field configurations are 
also considered, a classically stable ground state is possible at ipo 21 

Theories ( |2.1|) and ( |2.2| ) are conformally equivalent to general relativity with a 
matter sector containing two interacting scalar fields. 22 By redefining the graviton 
and dilaton fields as 

g^u = ft V> ^ = 2K 2 h(^) (2.5) 

and 



where k 2 
Hilbert form 



87rm P1 2 , the action (2.1) for U(ip) 



may be rewritten in the Einstein- 



S 



d 4 x 



A _ 1 ( v$) 2 - \AmVa) 2 - C{<S>)V{o) 



(2.7) 



where 

=n 2 = 2K 2 h(^P), C($) = A 2 ($) (2.8) 

and we assume that {h, dh/dip} > for consistency. The viscosity appears through 
the non-standard coupling, of the $-field to the inflaton's kinetic term. This 

coupling evolves towards unity as $ settles into a minimum of the effective poten- 
tial C($)V(cr), thus causing the accelerated expansion to end. In principle one can 
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therefore realise this scenario in general relativity if the matter sector of the the- 
ory is modified in an appropriate fashion. Some higher-dimensional, higher-order 
gravity theories exhibit this conformal structure upon compactification of the extra 
dimensions. 23 

In the following discussions we refer to g^ u as the Jordan-Brans-Dicke (JBD) 
frame and g^ u as the Einstein-Hilbert (EH) frame. In the former it is the matter 
contributions which are canonical, whereas gravity is canonical in the latter. The 
evolution of the Planck mass in the JBD frame is translated into a time-dependence 
for gauge and Yukawa couplings in the EH frame. 24 

2.2 Field Equations 



Extremizing the action ( |2.2[ ) with respect to arbitrary variations of the metric pro- 
duces the gravitational field equations 
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^[V^^-g^m], (2.9) 

where □ = g^ u V M V v and the energy-momentum tensor T^ u is defined as the func- 
tional derivative of the matter lagrangian including the inflaton field. 25 To proceed we 
assume a space-time with isotropic and homogeneous spatial sections (the Friedmann- 
Robertson- Walker universes) where the fields are functions of cosmic time only. In 
this case minimising ( |2.2p with respect to variations in a and \1/ leads to the equations 

dV 

a + 3Ha = — — (2.10) 
da 

and 

(V rr^\ dw^ 2 dU 

K = 2u[- + 3H-) +—— + —, 2.11 
\ ^ ^ / \ ^ / d^ ^ d^ 
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respectively. Combining Eq. ( 2.11Q with the trace of ( |2.9|) we arrive at the dilaton 
field equation 

(2lu + 3) (iff + 3H^) =2U- - + 8ttT, (2.12) 

v ' v ' d^ d# 

where T = ptotai ~ Sptotai is the trace of T^ v . Finally, the time-time component of 
gives the Friedmann equation 



3 (h> + JL) = t {% V + i£ _ + £ + I'M + ) . (2.13) 
V a 2 / 2 U 2^ * * V2 



where = —1, 0, +1 determines the spatial curvature. 

It is necessary to reheat the universe shortly after the expansion has become 
subluminal. The couplings of the inflaton to other relativistic matter fields should 
then become important because the potential is steep. However, reheating in these 
theories is difficult to model so we follow Morikawa & Sasaki 26 by introducing a 
phenomenological dissipation term into the field equations. It was shown that for an 
exponential inflaton potential and a Yukawa type coupling the inflaton field equation 
Q2.10Q becomes 27 

d^ 

a + 3H& = — C v a, (2.14) 

da 

where to zeroth order the dissipation 



d 2 V 

C v ~ fM a = /y— (2.15) 

is proportional to the effective mass of the inflaton and / = 0(1). This form for C v 
should be appropriate for any steep potential. The Bianchi identity V = then 
implies that 

p r = -AH Pr + C v a 2 , (2.16) 

where p r is the energy density of relativistic particles. The decay product must lead 
to baryogenesis, and we assume the reheat temperature is significantly larger than its 
rest mass mj. We may then treat it as a relativistic fluid in local thermodynamic 
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equilibrium. Once the inflaton's effective mass falls below mx the dissipation becomes 
negligible. 



2.3 Change of Independent Variable 

Numerical solutions must be found which trace the evolution of the universe from the 
Planck time to the current epoch. This corresponds to some 60 orders of magnitude 
in the independent variable. The codes used to produce these solutions become inef- 
ficient when integrating over more than 15 orders of magnitude. Therefore a change 
of independent variable was used to speed up the process of generating solutions. 

The number of e-foldings N = In a was used as the new independent variable since 
this is clearly a monotonically varying function for expanding universes. From the 
relation H = a/a = d(lna)/dt it can be seen that dt = dN/H. Hence for a variable 
x(t) we must make the substitutions: x — > Hx' and x — > H 2 x" + HH'x' where a prime 
denotes differentiation with respect to N. We now have the following equations for 
k = 0: 

„ / HH'\ 1 dV C v , . . 

a « + 3 + — U< + -- + ^/ = (2.17) 



p' r + p r - HC v a' 2 = 



(2.18) 



(2to + 3) 



+ 3 



n U * dU du , TA2 

2 1 1 (* ) 

H 2 H 2 d^ d^ 1 J 



•>7T 



AV(a) 

~1p~ 



a 



and 



dN 
If 



(2.19) 
(2.20) 



H 2 



1U 8vr /TA . 
2^ + T iV + p r> 



8tt a' 2 

+ 3 ¥ + W 



(2.21) 



We now need an extra equation to calculate H'. This is found by taking Q2.ll ) 
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and substituting for + 3H^>) from Q2.12 



) leaving 



2a; + 3/ * d^ 



+ 57T3?( 2C/ + to («'M-" 2 ))-(l) 2 - < 2 - 22 > 



or converting to independent variable N: 



6 V* 2cj + 37 V" d^ ^ + d^ 



+ 



1 _2o^_ 1 / 2[/ g?r Uy, ^ _ H 2 a 
6 2uj + 3 ^ V V v ; 



/2 



if 2 u; f ^) - 2H 2 . (2.23) 



Defining y = a' and z = we obtain a set of first order, ordinary differential 
equations which can be integrated given initial conditions on {t, p r , a, a', at 
some initial N — iVj. When certain slowroll conditions, such as ^ and 
(T 2 are valid, this coordinate transformation will simplify any analytical approach 
considerably 

Before proceeding to solve these equations for specific inflaton potentials, we shall 
investigate the most important observational constraints that any successful scenario 
of this type must satisfy We present a detailed discussion of how such limits arise 
since they will always be relevant in any future analysis. 

3 Constraints from the CMBR, nucleosynthesis and 
the solar system 
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3.1 Amplitudes of scalar and tensor perturbations in the con- 
formal frame 



The purpose of this section is to derive the general formulae for the amplitude of 



density fluctuations which arise in models based on theory (2.1) and then discuss the 
constraints which any scalar-tensor scenario based on a steep potential must satisfy. 
We generalize the calculation of Berkin and Maeda 28 by considering an arbitrary 
functional form for h{ip) or equivalently for oj(^>). If one wishes to employ Bardeen's 29 
analysis for the evolution of super-horizon sized perturbations, it is necessary for 
consistency to perform the calculation in the EH frame where the Planck mass is 
truly constant. To proceed analytically, however, it is also necessary to assume certain 
'slow-roll' approximations for the two fields and ignore all quadratic first-derivative 
and linear second-derivative terms in the field equations. The full equations derived 
from theory ( |2.7| ) were presented in Ref. 30 and for slowly rolling fields they reduce 
to 

H 2 « — CV (3.1) 
3 V ; 

3^$ » -C*V (3.2) 
2,Hd v a « -AV a , (3.3) 

where subscripts $ and a denote differentiation with respect to $ and o respectively, 
d^ = d/d?7, 7] = J dtQ(t) denotes cosmic time in the conformal picture and tildes refer 
to quantities defined in the EH frame. In general, the total energy density of the 
system is defined by 

p=\{^) 2 + \A{& n ay + A*V. (3.4) 

and the expression for the density spectrum can be found by extending the results 
of Lyth 31 to the case of two interacting scalar fields. If we denote by 771 and 7] 2 the 
times a perturbation leaves and re-enters the horizon in the EH frame, then 









p . 




p . 



(3.5) 



10 



where (3 = 2/[3(l + uj)} and uo = p/p is the ratio of the pressure to the energy density 
at the two epochs. In this discussion, the term 'horizon' refers to the inverse Hubble 
scale at the times rji and horizon crossing is defined in terms of comoving wavenumber, 
k(j]), by the expression k(r]) = a(i])H(i]). 
During inflation, u)\ ~ — 1 which implies 



~ S 1 
P 



1 12 



3 ll 



P + P 



(3.6) 



m 



It is necessary to derive an expression for 8p at the first horizon crossing. This is 
achieved by varying ( |3.4j ) and ignoring quadratic terms in d^$ and d v a. Dimensional 
analysis implies 8(d v $>) ~ H (£$) and S(d v a) m H5a, so 



5p « -2{Hd ri §5§ + A^HdrjaSa), 



(3.7) 



where the field equations (3.1) have been used to simplify the result. The terms 5$ 
and 8a are stochastic quantities arising from quantum fluctuations in the fields. In 
practise, one uses the two-point correlation functions to estimate their magnitudes. 32 
However, in this theory the inflaton has a non-standard kinetic term in the action 
fl2.7p . Naively, one would expect \8a\ ~ H, but an additional factor is present in this 
expression which arises when the second quantization is performed. In general, this 
factor is given by the inverse square root of the function coupled to the kinetic term 
of the inflaton. Hence 



By substituting 



\8®\&H, \8a\ w A- l/2 (<$>)H. 
into (13.61), we arrive at the final result 



(3.f 



A s = 


\5~p 


w aH 2 




_ P _ 


m 



[dggj + yC4|d^o-| 
(d (? $)2 + A(d v ay 



(3.9) 



where a is a numerical constant of order unity and the right-hand side is evaluated 
at the start of the last 60 e-folds of inflation. This is a very general expression and is 
valid for arbitrary functions {h(tp), V(a)}. 
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Two regions may be denned in terms of the relative values of d^ and y/Ad v a as 
Region I - \d„$\ < VA\d ri a\ (3.10) 



Region II - |d„$| > v^ld^cr! 



In these two limiting cases, the amplitude of density perturbations becomes 



(3.11) 



A s ,i 



A 3/2 



y3/2 

~v7 



A 



an 3 C 3 / 2 



5,11 



V l/2 



(3.12) 



(3.13) 



y/3 C$ 

where the field equations (3.1) have been used. In region II, the evolution of the 
dilaton field, as represented by the $-field, dominates the inflationary dynamics and 
one recovers the results for the original hyperextended scenario. 18 

The general formula for the perturbation spectrum in the EH frame can be ex- 
pressed in terms of quantities in the JBD frame with the use of equations ( |2.5| ) and 
Q2.6p . The conformal transformation fl2.5| ) implies that the two-point correlation func- 
tion of the inflaton in the JBD frame is \5a\ ~ H, as expected. It is straightforward 
to show that the expression for A s is given by 



a 



A l ' 2 H 2 



|$| + A^ld 



CM 



$ 2 + Aa 2 



(3.14) 



The results of McDonald 33 were derived in the JBD frame and they can be compared 
to those of Ref. 28 for the special case of the JBD theory by using Eq. ( |2.5| ). For the 
pure JBD theory, Eq. ( |3.14| ) reduces to 



Ac w aH 2 



(l + 6e)^ 2 + o- 2 



(3.15) 



where e = and this is identical to Eq. (47) in McDonald's paper with the 
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translation 

d/dt — > (1 + 6e) 1/2 d/dt (3.16) 

for the time derivative of the dilaton. Such a change is very close to unity when e C 1. 
Hence the change in Newton's constant is negligible compared to the evolution of the 
perturbations, as assumed by McDonald. (When e > 1/6, the slowroll conditions 
used to derive (3.1) break down and the derivation of Eq. ( |3~9l) is inconsistent). 

The close agreement between the expressions for the density spectrum in the 
two frames suggests that the conformal transformation (|2.5|) used to recast the JBD 
theory into the Einstein-Hilbert form may be valid at the semi-classical level in these 
chaotic models. The two results should be identical in the limit as e — > 0. A detailed 
comparison of the two different methods was made by Guth and Jain 34 within the 
context of old extended inflation in which the inflaton is fixed. These authors extended 
a previous analysis by Kolb et a/. 35 They also conclude that the technique of conformal 
transformations is valid up to a numerical factor of order unity when e < 1. 

Finally the expression for the amplitude of gravitational waves (tensor modes) 
should also be calculated in the EH frame. One can view a graviton as a minimally 
coupled massless scalar field with two degrees of freedom corresponding to the two 
polarization states of the wave. Abbott and Wise 36 have derived an expression for the 
amplitude at horizon crossing for an arbitrary inflationary solution, so their results 
are also valid for the theories under consideration in this work. Therefore, to a first 
approximation, the equivalent expression for the tensor modes is 

A o = ( 3 - 17 ) 

To summarize, we have derived the expressions for the scalar (As) and tensor 
(Aq) modes in the EH frame. For comparison with observation, however, we require 
the equivalent expressions in the JBD frame, since we are interpreting this as the 
physical frame. In general, the equivalent expressions As and Aq are related to their 
counterparts in the EH frame by an expression involving the conformal transformation 
( |2.5|) . The relation therefore depends on the form of h(tf)). However, if the condition 

d„Q 

< 1 3.18 
MI V ; 
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holds, and if fi is a function of t only, it is straightforward to show that the scale 
factors and expansion rates in the two universes are related by 

a(t) = JT^afa), H = QH. (3.19) 

Eq. (|3.19|) implies that k = aH « aH = k, so the scale dependences in the two frames 
are approximately equal. Moreover, the definition of the energy density, p = g 00 T 00 , 
implies 

55 5p SQ , 
JL = JL-±— 3.20 
p p si 

and we conclude that amplitudes are also equal when Sfl/Q -C 1, as implied by Eq. 

for drj fa H~ 137 Therefore, it is sufficient to consider results from the conformal 
(EH) frame directly. 

3.2 The CMBR and large-scale structure 

If particle physics specified the unique inflaton potential, the amplitude A s would be 
known for all scales. Because there exist many possible models, however, the ampli- 
tude of fluctuations must be normalized using observations of large-scale structure. 
The CfA survey suggests that the rms fluctuation in the galaxy counts is unity in 
a sphere of radius 8/i _1 Mpc and we normalize the amplitude by specifying the rms 
fluctuation in the mass distribution to be cr 8 = b% 1 within this sphere. 53 6§ is the 
biasing factor which we assume to be constant over the scales of interest. 

We shall show in section (5) that the inflationary solutions approximately 60 e- 
foldings before reheating can be expressed as a power law a(t) oc t p for p ^> 1. The 
slowroll approximations are valid in this region of parameter space. It is well known 
that this solution leads to a primordial power law spectrum of scalar fluctuations, 
P(k) oc A 2 s (k)k oc k n , where n — 1 — 2/(p — 1) is the spectral index. The spectrum 
of tensor modes has an identical scale-dependence for this solution when n < 1, i.e. 
Aq oc As, and scale-invariant fluctuations correspond to n — 1 (p — oo). 38 

Equations relating <7g to n have been derived using results from the CMBR and 
large-scale structure. 12 The range of parameter space (crg,n) consistent with these 
observations can then be determined. 39 For our purposes, these constraints restrict 
the value of the effective JBD parameter ( |2.3[ ) during inflation. 
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At present the observation of multipole anisotropies by the COBE DMR exper- 
iment provides the strongest constraint from the CMBR. The root of the variance 
at 10° is observed to be a T (10°) = (1.1 ± 0.2) x 1(T 5 , where the mean blackbody 
temperature T = 2.736 K is taken and limits are to 1-sigma. 12 On 10° scales reion- 
ization processes such as Thompson scattering are not important and the dominant 
contribution to the anisotropy is from the Sachs- Wolfe effect when photons at decou- 
pling are redshifted as they climb out of potential wells. 41 For power law inflation the 
observed anisotropy is due to both scalar and tensor modes. 42 Following an identical 
procedure to Ref. 40, a numerical fit relating a^(10°) to n and cr 8 can be found and 
equated to the observed 10° variance to yield 

a 8 = (1.18 ± 0.2)e 2 - 63 ( 1 - n )^(3 - n)/(15 - 13n). (3.21) 

The contribution of the tensor modes arises solely in the rooted factor and it becomes 
negligible as n approaches unity. It has the effect of decreasing a 8 for a given spectral 
index, thereby increasing the allowed bias. 

A second result was presented in Ref. 39 based on the IRAS/QDOT 13 and 
POTENT 43 galaxy surveys. Quoting a value of 6/ = 1.2 ± 0.3 at the 2-sigma level for 
the IRAS biasing factor, these authors deduced that 



19 \ 2 / 3 

a 8 = (0.91 ± 0.25) ( '■ . (3.22) 



.2.9 -n, 

Eqs. ( |3.21| ) and ( p.22|) can be combined to yield a bias-independent lower limit on the 
spectral index of n > 0.84 for consistency between COBE and QDOT. It is this limit 
which most strongly constrains the JBD parameter during inflation. If Eq. (3.11) is 
valid in the pure JBD theory one recovers the old extended inflationary solution. 44 
In this regime the spectra have the power law form discussed above with a spectral 
index 

n = l -4 • (3.23) 

Consistency with the CMBR and QDOT results therefore requires u > 25 at the start 
of the last 60 e-foldings before the end of inflation. It is important to note that the 
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constraint due to bubble collisions, u < 18, does not apply here because the phase 
transition is second-order. 19 



3.3 Nucleosynthesis and solar system constraints 

After the expansion becomes subluminal the most important constraints on the value 
of uj and the vacuum energy arise at the nucleosynthesis era and the current epoch 
respectively. The standard model of primordial nucleosynthesis can be used to con- 
strain the vacuum energy density and the JBD parameter at temperatures 0(1) MeV. 
Time-varying gravitational and cosmological constants modify the expansion rate of 
the universe and therefore the nuclear reaction rates during this epoch. At tempera- 
tures exceeding the 'freeze-out' temperature Tp m 0.8 MeV the neutrons and protons 
are held in local thermodynamic equilibrium by weak interactions. 15 In the standard 
big bang picture the neutron-to-proton ratio at T > T F is determined by the Boltz- 
mann factor (n/p) = exp(— Q/T), where Q = 1.293 MeV is the n-p mass difference. 
The energy density of the universe is then p = n 2 g e sT 4 /30, where g e s is the number of 
relativistic degrees of freedom. The nucleons fall out of equilibrium when the reaction 
rate equals the expansion rate of the universe, i.e. 



The temperature drops below Tp and free neutrons undergo /9-decay until the pho- 
todissociation of deuterium becomes energetically unfavourable at some temperature 
T D . At this point synthesis of He A proceeds rapidly. 

The abundance of ife 4 depends crucially on the number of neutrons present at 
Tp,. There will always be two competing effects in any modification to the standard 
picture. Increasing the freeze-out temperature increases (n/p) and one might expect 
the He A abundance to be correspondingly higher. However, in many cases it then 
takes longer for the universe to cool to T — T D . More neutrons can undergo /5-decay 
and this effectively reduces the He 4 abundance. 

When vacuum energy is introduced the second effect is dominant and the observed 
He A abundance therefore leads to an upper limit on VLy. lA From detailed numerical 




(3.24) 
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calculations it was found that this limit is 

n v [lMeV] < 0.08 (3.25) 

for three neutrino species. 14 

There is also a limit on the change in Newton's constant with time. In the standard 
JBD theory the strength of gravity is higher at earlier times and this increases Tp as 
indicated by Eq. ( |3.24| ). One therefore expects an upper limit on G/G to exist or 
equivalently a lower limit on u. Casas et al. 45 found u > 250, but by dropping some 
of the simplifying assumptions used to obtain this bound, Serna et al. 46 found the 
weaker limit of uo > 50. As a first approximation we shall apply this weaker limit to 
the more general scalar-tensor theories under consideration here, i.e. 

lMeV] > 50. (3.26) 

It should be emphasized that this limit is only strictly valid when uo is constant for all 
time. Although constraints ( p. 25] ) and ( 3.26 ) may be weaker than those derivable from 



more exact calculations, they are sufficient to severely limit the scenario discussed 
here, as is shown by the numerical calculations in section 4. 



Finally, constraints on the magnitude and form of the JBD parameter (23) at 
the present epoch can be obtained from solar system experiments by using the post- 
Newtonian approximation. In this analysis one considers the time independent spher- 
ically symmetric metric around a point mass m, expanded as a series in the gravita- 
tional potential U = m/r, i.e. 

ds 2 = (1 - 2aU + 2/3U + ...)dt 2 + (1 + 2 7 £f)dx 2 , (3.27) 

where dx 2 = dxidx 1 (i = 1, 2, 3). Current observational limits on the Post-Newtonian 
parameters are 

a = l±10 -4 , (7 + l)/2 = 1 ± 10~ 3 , (2 + 2 7 - 0)/3 = 1 ± 10~ 2 , (3.28) 

whereas general relativity predicts the values a = (3 = 7 = 1. Nordtvedt 16 analyzed 
the generalized scalar-tensor theory (|2.2|), finding that 



'- i+ ( « + SgW (3 - 29) 
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and 

7 = ^^ =>- u#o,3K]> 500. (3.30) 

It is important to note that uj — > +00 is a necessary but not sufficient condition for the 
recovery of general relativity at the present epoch. We also require from Eq. ( |3.29| ) 



that (do; / d$f) / 'uj 3 — > as uj diverges. However, this condition is always satisfied when 
the dilaton is located within the vicinity of a turning point in h(ip), because 



new 



To summarize, the most important limits are uj > 25 during inflation, uj > 50 at 
nucleosynthesis and uj > 500 at the present epoch. 

4 Numerical Results 
4.1 Initial Conditions 

A number of plausible models were considered by numerically integrating the system 
of equations ( |2.17| )-( [2T2~1"| ) plus ( |2.23| ). The analysis does not depend strongly on the 



precise form of h(ip). One only requires that it be at least C 2 continuous, contain a 
turning point and be normalized to unity at this point. In this sense the scenario we 
are investigating is rather generic. To ease the calculation we specified 

\Er = h(if)) = sin 2 Pip, (3 = constant, (4.1) 

because this leads to a simple form for ou(ty) given by 

(jfijf) = — _ — . (4.2) 

Although there are no known particle physics models which lead directly to Eq. ( |4.1|) , 
it has a simple analytical form which can be viewed as an approximation to a more 
complete theory. Indeed, for small ^ one can treat Eq. ( |4.2|) as a perturbation 
to the JBD theory up to terms including 0{^ 2 ). General relativity is recovered 
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as \I/ — > 1. Following Linde, the initial conditions were specified by treating the 
quantum boundary (QB) as the most natural set of initial conditions for chaotic 
inflation in scalar-tensor theories. 47 This boundary corresponds to the hypersurface 
at the Planck density and represents the earliest times at which initial conditions 
can be placed on classical fields. In theory ( |2.1[ ) the Planck mass is related to the 
dilaton by m^f) = h(ip) so the QB is reached when V((Ji) ~ h 2 (ipi) = mpj. There 
is therefore a 1-jet family of initial conditions relating the values of a and ip. Indeed, 
the probability that the universe can be created from nothing via quantum tunneling 
is 

p^exp(-3m 4 1 (^)/81/(or)), (4.3) 

which is only high on the QB. 48 

In the numerical calculations we employ this argument as a first approximation. To 
be more accurate, though, one should also consider the effects of quantum fluctuations 
in the two fields. Once it has started inflation always occurs at some point in the 
universe if the change in the fields due to quantum fluctuations with wavelength larger 
than H~ x exceeds the classical change in the time interval H~ x due to the equations of 
motion. 47 An island universe resembling our own (almost) flat Friedmann space-time 
may only form when such a situation is reversed, and the point of equality is a more 
accurate estimate of initial conditions. 

The evolution of the quantity dlna/dlnt = Ht with respect to the number of 
e-foldings, N = In a, was investigated because this allows a number of interesting 
features to be shown diagrammatically. For example, if the scale factor expands as a 
power law, a oc t p , Ht = p = constant and the graph is a horizontal line. On the other 
hand, the exponential (de Sitter) expansion may be written as N = Ht + constant 
and this is a straight line of gradient 7r/4. The critical solution for inflation is the 
Milne universe, Ht = 1, and the strong energy condition is violated above this line. 
The end of inflation can therefore be defined as the point where the graph cuts the 
line Ht = l. 
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4.2 Decaying Power Law Potentials 

Witten has shown that potentials of the form V oc <J~ a , where a > 0, can arise when 
super symmetry is spontaneously broken. 49 The symmetry breaking occurs when the 
potential has non-zero values. Therefore, one would have a natural candidate for 
the inflaton if a was the field responsible for breaking supersymmetry. In Einstein 
gravity the potential is too steep at small a to lead to inflation, but at large a an 
intermediate inflationary solution a oc exp(Ft q ) is found, where F and q < 1 are 
positive constants. 50 

The question we address is whether inflation at small a is possible in the scalar- 
tensor theory discussed above and the numerical results are shown in Figure (1) for 
a = 10 and (3 = 0.07. A power law expansion with p = l/(8(3 2 ) + 1/2 arises when 
100 < N < 260. In this region u(^f) is approximately constant (\1/ 2 <C 1) and 
the solution corresponds to that found in the JBD theory when the inflaton is held 
constant. This is the old extended inflationary solution a oc t^+V 2 . 44 Hence the known 
results for the form of the perturbation spectra derived when Eq. (3.11) is valid can 
be employed in this analysis and Eq. ( P-13j ) applies. Consistency with the CMBR 
and QDOT results therefore requires (3 < 0.07. 

Figures (la) & (lb) 

Unfortunately there is a graceful exit problem in this model. Figures (la) and 
(lb) show how the solution temporarily approaches the Milne limit as a 2 increases 
and Einstein gravity is recovered at N « 310. At this point the intermediate solution 
takes over where Ht = qN + constant. The gradient of this line is consistent with the 
analytical result a = 4(g _1 — 1) derived by Barrow. 50 (This provides a useful check for 
the numerical calculations.) Because the de Sitter solution is an attractor at infinity 
in this case, it is not possible to obtain subluminal expansion unless the reheating 
process is highly efficient. Rather than redshifting to zero, the radiation density grows 
during the power law phase due to the continual decay of a. This is a generic feature of 
these models, but is not sufficient to establish a radiation-dominated phase as required 
by nucleosynthesis. Naively one may think that & 2 would increase sufficiently fast as 
^ —+ 1, but it is too small relative to ^ during inflation to contribute significantly to 
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the reheating. We considered models up to a = 100 and found the same qualitative 
behaviour. 



4.3 Exponential potentials without dissipation 

We proceeded to investigate the exponential potential V(a) = V exp(—\K<j) where 
{Vo, A} are independent constants and k 2 = 8tt. In Einstein gravity this potential 
leads to an attractor power law solution a oc t 2 ^ 2 when A 2 < 6 and inflation occurs if 
< A 2 < 2. The attractor solution becomes a oc t 1 ^ 3 for all A 2 > 6 and the physical 
reason for this is as follows. In the spatially flat Friedmann model the A 2 < 6 attractor 
arises because the kinetic and potential energies of a redshift at the same rate. This 
is the deep reason why inflation can never end in this model. When A 2 > 6, however, 
the potential is too steep and the field becomes effectively massless as V — > 0. 51 

A number of new features arise when the Planck mass is allowed to vary. We find 
that inflation is possible when A 2 > 6 and that the expansion becomes subluminal 
without any fine-tuning. As an example, we discuss results for A 2 = 12 and /3 2 = 1/8. 
In Figure (2) the dissipation term Cy is removed for all time. The p = to + 1/2 
power law solution is again found, thereby implying that the dilaton dominates the 
dynamics. The expansion does indeed become subluminal and settles around the 
value p = 1/2. This result is independent of the value of j3 and is understood by 
investigating Eq. (|2.12| ). For U = and u approximately constant, this equation can 
be rewritten in the form 

q + 3HV + —^^0, (4.4) 
where V e g is an effective interaction potential defined by the integrability condition 

= T (4.5) 

2cj(^) + 3 v ' 



and T = 4V — a 2 . Eq. ( |4.4j ) resembles the equation of motion for a minimally coupled 
field and we may view \1/ as evolving along its effective interaction potential. During 
the inflationary epoch, V ^> a 2 and has negative gradient. Hence the values of 
\1/ and a 2 increase. This means that T, and therefore dV e ff/dty, eventually change 
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sign as 4V — & 2 passes through zero, thus causing the dilaton to slow and reverse its 
motion. Hence, there exists a damped oscillatory behaviour around 4V = a 2 and some 
constant value of The theory soon resembles a rescaled version of general relativity 
with a traceless energy-momentum tensor; the solution is therefore equivalent to the 
radiation-dominated universe a oc 



Hence, the effect of a varying Planck mass is to change the attractor for steep 
exponential potentials from p — 1/3 to p — 1/2. The above argument should apply 
to any steep potential for which T < in Einstein gravity. Including an uncou- 
pled radiation component will not alter the argument since its energy-momentum is 
identically traceless and does not affect V c g. 

4.4 Exponential potentials with dissipation 

In Figure (3) the dissipation term is included for all time. Figure (3c) shows the 
radiation density growing during inflation and a sufficient reheating temperature for 
baryogenesis to proceed is easy to obtain in this model. Moreover, from Figures (3b) 
and (3d), it is seen that the radiation begins to dominate once Einstein gravity is 
recovered and the dissipation prevents the quantity 4V — a 2 from passing through 
zero. However, a new scaling solution is found, where p a / p ra ,A = constant, as is shown 
in the late time behaviour of Figures (3a) and (3d). 



This solution may be derived analytically for the special case where Cy is constant. 
We consider a model in which a particle species X, with equation of state px = 
(7 — l)px for some constant 7, is coupled to the inflaton through an interaction 
p x = Tx<7 2 , where Tx is constant. The Bianchi identity and Friedmann equation for 
this system are 



where py represents the inflaton energy density and we assume Einstein gravity holds. 



Figures (2a) & (2b) 



Figures (3a) - (3d) 



(1 + 2T x )t + 3(1 + -fF x )H& + V a = 
3H 2 = k 2 ( Pv + px) = k 2 + Tx) & 2 + V , 



(4.6) 
(4.7) 



22 



By differentiating Eq. (|4.7[) with respect to cosmic time and substituting Eq. 
for a we arrive at the 'momentum' equation 



K 



^ = -y(l + 7rx)a 2 . 



(4.8) 



This implies that Eq. ( [4.7] ) may be rewritten in the Hamilton- Jacobi form by defining 
a new time coordinate 

t = ~(l -I'v) / da'f^rl . (4.9) 



: 1+ ,r,)/^g§r 



We find 



3n 2 (l + ^xfH 2 - 2(1 + 2Y X ) (^pj = k\1 + 7 r x ) V, 



(4.10) 



which yields the attractor solution 



H = \/]4exp (-Akct/2) , A/«x = 21n 



21 + 7^ 



(4.11) 



where v4 is a positive constant. Hence the scale factor grows as a power law a oc t p , 
where 

p = A(i + 7 r x ). (4.12) 

(When Fx -> we recover the standard power law solution with p = 2/A 2 ). The 
contribution of X to the matter content of this universe may be expressed through 
the quantity Q x = K 2 px/3H 2 = px/(px + Pv)- F° r the solution (|4.11|) we find 



ft 



Tx\ 2 



x 



3(i + 7 ry 



(4.13) 



is constant, implying Eq. (|4.11| ) represents a scaling solution and the matter and 
radiation densities redshift as py oc px oc t~ 2 oc a~ 2 / p . Moreover, by substituting Eq. 
( [4.1 1| ) into Eq. ( [4.13|) , the condition Vq > implies this attractor exists only if 



2(1 -ft 



(4.14) 



x, 
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The nucleosynthesis constraint ( |3.25|) then implies Tx > 23/4 which yields the lower 
limit A 2 > 36. 

In principle this result suggests the nucleosynthesis constraints can be satisfied in 
this model if A 2 is sufficiently large. However, if the dissipation is not removed, the 
scaling solution will survive through to the decoupling era. Numerical results have 
shown that the lack of observed spectral distortions in the CMBR imply Qy < 4x 10~ 4 
if the vacuum decays into low energy photons. 14 For consistency, Eq. ( 4.14 ) then 
implies A 2 > O(10 3 ), which is clearly unrealistic. In any case it is physically reasonable 
to suppose the dissipation becomes negligible once the effective inflaton mass falls 
below the rest mass of its decay product. Naively one would expect the vacuum 
energy to rapidly redshift to zero at energies below this mass scale, with the radiation 
density falling as a -4 . In general relativity this would be the case, but the evolution 
of the Planck mass again significantly alters the arguments, as shown in Figure (4). 

Figures (4a) & (4b) 

As soon as the dissipation term was removed the vacuum rapidly dominated the 
universe once more. This feature arises because the peak of reheating occurs while 
the dynamics is still dominated by the dilaton's motion. The dilaton viscosity still 
slows the inflaton and thus its energy density does not decay as fast as it would in 
Einstein gravity. On the other hand, the time-dependence of the Planck mass does 
not affect the radiation, which still decays as p r oc a -4 . This follows because the 
energy-momentum tensor of the free radiation field is traceless and does not appear 
in the dilaton field equation (2.12). We suggest that this qualitative behaviour should 
not depend too strongly on the precise form of the potential unless V(a) is very steep. 
In this case, though, the vacuum energy would be insignificant at the present epoch. 
The same qualitative behaviour was observed for values of uj = 500 (/3 = 0.0158). 
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5 Conclusions and Implications 



The philosophy behind this work was to identify a decaying cosmological constant 
at the present epoch directly with the vacuum energy which drove the inflationary 
expansion without altering the form of the potential. This would solve, without 
severe fine-tuning, the cosmological constant problem and could provide a possible 
explanation for a number of observational results. The limiting solution for acceler- 
ated expansion, a oc t, arises when the vacuum energy density varies as pv oc a~ 2 . 
Hence, pv must redshift faster than a~ 2 if inflation is to end, but must decay slower 
than the pressure-free matter component, p m atter oc a -3 , if it is to dominate at late 
times. This requires viscosity to be present at early times. We derived expressions 
for the amplitudes of scalar and tensor fluctuations with a general inflaton poten- 
tial and scalar-tensor theory, and the strongest constraints on any model arise from 
the CMBR, primordial nucleosynthesis and solar system observations. These con- 
straints are important in any scenario of this type. In general, it is difficult (if not 
impossible) to satisfy these constraints simultaneously if the viscosity arises due to a 
time-dependent Planck mass. 

We considered two simple forms for a steep potential, V oc a~ a and V oc exp(— \<j). 
In the former a graceful exit problem arises for realistic values of a. The latter is 
more promising and a number of scaling solutions were found both numerically and 
analytically. In particular, inflation occurs and ends naturally when A 2 > 2. This 
model leads to a power spectrum consistent with observation, but the vacuum does 
not decay sufficiently fast to satisfy the nucleosynthesis constraint ( pj.25| ). It also ap- 
pears that this constraint cannot be satisfied for very steep potentials, which decay 
rapidly to zero in Einstein gravity, because the Planck mass is still evolving during 
reheating. This causes pv to redshift at a much reduced rate relative to the radiation 
component. 

To summarize, the scenario as outlined in section 2 is not viable for the examples 
considered here and requires more complicated potentials and further extensions. 
However, with this numerical code, it is possible to develop working hyperextended 
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chaotic models based on the theories discussed in section 2 if the steep potential has 
a minimum. One extension to this analysis is to investigate any effects a dilaton self- 
interaction potential U(ip) may have. Alternatively one could consider more general 
couplings of the dilaton and inflaton fields in action ( |2.1| ) or alter the form of the 
dissipation term Cy which models the reheating process. 

These results may have implications for the joint evolution of the comological and 
gravitational constants. A present-day vacuum term may arise if 

1. The inflaton settles into a minimum of its potential at V = and a second 
scalar field is located in a non-zero minimum of its own potential. 

2. The global minimum of the inflaton potential is located at V ^ 0. 

3. V(a) decays monotonically for all time. 

The third possibility is the most attractive, but based on the above numerical 
calculations, one may conjecture that the vacuum will generally not decay sufficiently 
fast relative to the relativistic matter if rn.pi is time-dependent during and shortly 
after inflation. Our analysis therefore favours a reheating process via oscillations of 
the inflaton about some minimum if the Planck mass varies with time, but it does 
not rule out other possibilities. This agrees indirectly, and for different reasons, with 
the conclusions of Ref. 52, in which a detailed quantum mechanical description of 
the reheating process was given. These authors concluded that particle production is 
only significant during the oscillating phase of the inflaton. (A direct comparison of 
conclusions cannot be made, however, since these authors only considered reheating 
in Einstein gravity, whereas the dilaton is still evolving during reheating in the models 
discussed here.) 

Conversely, if one prefers a continually decaying potential, this suggests a constant 
Planck mass is required. In any case, it appears that some degree of fine-tuning is 
necessary if a cosmological constant arising from the inflaton potential is to be non- 
zero at the present epoch. 
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Figure Captions 

Figure 1: Numerical solutions for the potential V(a) = 10~ 1( V~ 10 , with (3 = 0.07 
and initial conditions a = ^ = 0. The dissipation coefficient / = 1.0 and the rest 
mass of the decay product X is m x = 10 9 GeV. Figure (la) illustrates the evolution 
of the power index of the solution. Power law and intermediate inflationary expansion 
is observed, but a graceful exit from the inflationary regime is not found. Figure (lb) 
illustrates the evolution of the a;-parameter in terms of tan _1 u;(\I/). The intermediate 
solution takes over once u> diverges to infinity (i.e. tan _1 c<j — > +7r/2) and Einstein 
gravity is recovered. 

Figure 2: Numerical solutions for the potential V oc exp(— \kg) with A 2 = 12, 
/3 — 1/8 and initial conditions a = ^ = 0. Dissipation effects have been removed 
and / = 0. Figure (2a) shows that inflation is possible and ends naturally as the 
expansion becomes subluminal. The expansion approaches the attractor p — 1/2 
as u(ip) settles to the constant value in Figure (2b). A rescaled version of general 
relativity is recovered but the Planck mass is smaller than the currently observed 
value. 

Figure 3: The same conditions as in Figure (2), but with f3 = 0.025 and / = 1.0. 
The vertical dashed line measures where the expansion becomes subluminal. Figure 
(3c) plots the evolution of the radiation density (p r ) and the inflaton energy density 
(p a ). Figure (3d) plots the ratio of the inflaton energy density to the total energy 
density. The effects of dissipation are included for all time. A new scaling solution 
is found once Einstein gravity is recovered, but this violates spectral distortion limits 
on the CMBR. 

Figure 4-' As Figure (3), but the dissipation terms are removed at m x = 10 10 GeV. 
In Figure (4a) the inflaton rapidly dominates the expansion after the dissipation is 
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removed. This follows because the dilaton is evolving during the reheating, as shown 
in Figure (4b) by the evolution of u(^f). 
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This figure "figl-l.png" is available in "png" format from: 
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